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ANALYTICAL  STUDIES  PERTAINING  TO  THE  MEASUREMENT  OF  LIGHT 


ABSORPTION  IN  AIR 

Bela  A.  Lengyel 

San  Fernando  Valley  State  College 


ABSTRACT 

Experimental  methods  proposed  for  the  measurement  of  the  absorption 
coefficient  of  light  in  air  are  analyzed.  The  calculations  are  primarily 
intended  for  the  0.5  to  1.2  micron  wavelength  region  where  for  certain 
wavelengths  the  absorption  coefficient  to  be  measured  is  estimated  to  be 
of  the  order  of  lO-^  cm-1.  (For  the  10.6  v  wavelength,  the  absorption 
is  Icrger  and  the  scattering  is  less  of  a  problem.)  The  effects  of  scat¬ 
tering,  heat  conduction,  static  and  acoustic  pressure  variations  are  ex¬ 
amined  and  are  related  to  the  problem  of  detecting  the  signals  produced 
by  the  absorption  of  heat  from  the  laser  beam.  Both  pulsed  and  cw 
schemes  are  considered.  The  measurement  of  the  static  pressure  rise  in 
a  narrow  tube  appears  as  a  feasible,  although  marginal  method  which 
offers  some  hope  of  success . 
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Analytical  Studies  Pertaining  to  the  Keasurement  of  Light  Absorption  in  Air 


1.  Statement  of  the  Purpose  of  the  Report. 

It  is  of  interest  to  know  the  amount  of  heat  evolved  in  a  material 
through  which  a  laser  beam  is  passing.  A  parallel  beam  of  light  is  attenu¬ 
ated  as  it  passes  through  ma  .ter  and  this  attenuation  is  mainly  attributed 
to  two  effects:  scattering  and  absorption.  Scattering  on  material  parti¬ 
cles  results  in  the  directional  redistribution  of  the  energy  of  the  incident 
electromagnetic  wave.  Absorption  results  in  the  conversion  of  electromag¬ 
netic  energy  into  other  forms  of  energy.  The  energy  lost  by  the  electro¬ 
magnetic  field  is  ultimately  converted  into  heat.  This  means  that,  possibly 
through  the  intervention  of  processes  of  atomic  and  molecular  excitation, 
the  energy  abstracted  from  the  electromagnetic  field  is  redistributed  among 
the  available  degrees  of  freedom  of  the  material  (gas),  so  that  it  is  ulti¬ 
mately  undistinguis'nable  from  energy  obtainable  from  contact  with  a  heat 
rsservoir.  We  shall  assume  that  the  atomic  and  molecular  processes  res¬ 
ponsible  for  the  redistribution  of  energy  among  the  available  degrees  of 
freedom  take  place  so  fast  that  at  the  scale  at  which  we  contemplate  our 
observations  the  conversion  of  the  absorbed  energy  into  heat  may  be  regarded 
as  instantaneous.  It  should  be  noted  that  we  are  talking  here  about  energy 
redistribution  locally  and  not  about  its  redistribution  in  space.  The  latter 
generally  proceeds  at  a  much  slower  rate.  We  shall  thus  proceed  under  the 
assumption  that  the  energy  absorbed  from  the  electromagnetic  field  instantly 
appears  at  the  site  of  absorption  in  the  form  of  heat.  Such  an  assumption 
is  probably  valid  if  our  time  resolution  is  no  finer  than  10- 7  sec. 
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The  absorption  of  light  in  solids  and  liquids  is  usually  significant 
enough  so  that  a  direct  measurement  of  the  absorption  is  not  difficult.  In 
gases,  such  as  pure  air,  the  measurement  of  the  absorption  of  .light  is  very 
difficult,  because  the  rate  of  absorption  is  small  and  in  the  visible  region 
it  is  completely  masked  by  the  much  larger  effect  of  scattering.  The  situ¬ 
ation  is  somewhat  more  favorable  in  the  infrared  beyond  1.2  microns. 


The  optical  effects  of  scattering  and  absorction  combine  so  that  he 
intensify  of  a  plane  wave  passing  through  matter  is  described  by  the  equation 

I(x)  =  Ioe-S*  ,  (1.1) 


where  the  attenuation  coefficient  8  is  the  sum  of  two  terms: 

8  =  a  +  a.  (1.2) 

Here  a  is  the  absorption  coefficient,  o  the  scattering  coefficient  and  x  is 
the  direction  of  propagation. 


In  the  spectral  region  where  the  ratio  a/a  is  large,  it  is  not  feasible 
to  determine  a  by  measuring  the  attenuation  of  a  light  beam  through  matter. 
The  value  off  must  thus  be  determined  not  by  the  change  Droduced  in  the  ra¬ 
diation  but  by  the  change  produced  in  the  absorbing  matter.  We  already  noted 
that  the  absorption  of  radiation  results  in  an  increase  in  the  thermal  energy 
°f  the  absorbing  matter,  therefore  we  shall  address  ourselves  to  problems 
pertinent  to  the  detection  and  measurement  of  small  quantities  of  thermal 
energy  deposited  in  the  path  of  a  narrow  light  bqam  of  the  type  that  may  be 
obtained  from  a  laser.  The  spectral  region  where  a/a  <  1  is  not  included 
in  our  discussion  because  scattering  is  less  of  a  problem  and  the  expected 
values  of  a  are  greater  for  the  interesting  wavelengths. 
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It  has  been  estimated  that  for  certain  wavelengths  in  and  arouna  the 
visible  region  of  the  spectrum  a  is  of  the  order  of  10“9  cm-1.  Scattering 
in  pure  air,  as  calculated  from  Rayleigh's  theory,  produces  a  much  higher 
attenuation.  For  green  light  (0.5  microns),  a  =  1.7  x  10~7  cm  for  air  free 
of  all  contaminants.  We  reserve  the  quant. itativi  data  for  the  next  section, 
but  we  note  here,  tnat  around  10  microns  the  situation  is  quite  different. 
First,  significant  and  spectrally  selective  absorption  is  likely  to  he  en¬ 
countered  because  of  the  presence  of  the  molecular  bands  of  the  constituents 
cf  air  in  this  region.  Second,  Rayleigh  scattering  on  molecules,  which 
varies  with  the  fourth  power  of  the  wavelength,  is  reduced  to  negligible 
values  around  10  microns. 

For  the  reasons  explained,  it  is  best  to  regard  these  studies  as  pri¬ 
marily  applicable  in  the  spectral  region  0.5  to  1.2  microns  and  intended  as 
a  background  for  experiments  contemplated  with  ruby,  helium-neon  and  neo¬ 
dymium  las°rs.  Although  many  of  the  conclusions  reached  are  valid  over  a 
much  broader  spectral  range,  it  is  unlikely  that  they  are  relevant  insofar 
as  they  pertain  to  practical  methods  of  measuring  o.  For  the  10.6  y  wave¬ 
length  of  COj  lasers,  it  is  expected  that  a  is  greater  than  10~9  cm-1  — 
hence  should  be  easier  to  measure — especially  since  Rayleigh  scattering  s 
negligible. 

This  report  contains  calculation  pertaining  to  thermal  and  acoustic 
effects  expected  as  a  result  of  the  passage  of  a  powerful  laser  beam  through 
air.  The  coni igurations  examined  are  such  which  are  likely  to  be  chosen  in 
experiments  aimed  at  the  determination  of  a  in  the  0.5  to  1.2  micron  spectral 
region.  We  have  included  here  the  analysis  of  the  pressure  change  in  a 


completely  enclosed  cylinder  whose  wells  ere  kept  at  a  constant  temperature, 
the  calculation  of  the  rate  of  dissipation  of  heat  arc.  1  a  laser  beam  by 
means  of  thermal  conduction,  the  generation  of  a  quasi- static  pressure  in¬ 
crement,  and  an  acoustic  pulse  by  a  laser  pulse  and  some  problems  pertaining 
to  the  detection  of  an  acoustic  pulse  or  a  sequence  of  acoustic  pulses.  The 
material  as  acsembled  here  does  at  represent  a  complete  investigation.  It 
is  a  group  of  exploratory  calculations  preliminary  to  the  experimental  de- 
termination  of  the  absorption  coefficient  of  light  in  air. 


2.  Relevant  Properties  of  Air. 

IR  tHlS  section  we  have  collected  such  physical  data  concerning  air 
that  seem  relevant  for  computations  in  later  sections,  or  that  may  be  needed 

in  evaluating  an  experimental  scheme  proposed  for  the  measurement  of  the  ab¬ 
sorption  coefficient  a. 

Atmospheric  air  free  of  dust  and  other  contaminations  is  a  mixture  of 
gases,  primarily  0,,  Ar,  COj  mid  water  vapor.  Other  gases  are  present 
in  such  small  proportions  that  their  presence  does  not  noticeably  affect  the 
physical  properties  of  air  that  are  of  concern  to  us.  Of  the  gases  listed 
here  only  the  proportion  of  water  vapor  is  variable. 

In  what  follows  we  shall  assume  that  we  are  dealing  with  clean,  dry 
atmospheric  air.  By  this  substance  we  mean  a  gas  mixture  from  which  all 
water  vapor  has  been  removed  and  which  contains  the  gases  N2,  02,  Ar  and 
C02  in  the  i  roportion  as  they  are  universally  found  in  the  atmosphere1. 

1.  Handbook  of  Chemistry  and  Physics,  46th  ed.  p.  F-116. 
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We  note  in  particular  that  such  air  contains  0.033  percent  CO  by  volume. 
Although  this  proportion  seem;,  small,  it  is  we)l  to  remember  that  a  rela¬ 
tively  small  amount  of  CO  or  water  vapor  significantly  affects  the  op¬ 
tical  properties  of  air. 

It  is  customary  to  specify  the  physical  properties  of  gases  at  stan¬ 
dard  temperature  and  pressure;  i.e.  at  0°C  and  at  760  torr.  At  these 
conditions  the  properties  of  dry  atmospheric  air  relevant  to  this  study 
are  as  follows: 

Density:  pQ  =  1.293  x  10- 3g  cm-3. 

Specific  heats:  cp  =  0.240,  cv  =  0.172  cal  g-1  °C-1, 

Sound  velocity:  Co  =  3.31  x  104  cm  sec-1. 

Heat  conductivity:  K  =  5.68  x  10-5  cal  cm-1  sec-1  0C_1, 

Index  of  refraction  for  X  =  5000  A:  n  =  1.000294. 

At  standard  temperature  and  pressure  the  particle  density  of  every 
ideal  gas  is  L  =  2.69  x  10 19  cm-3.  (Loschmid's  number.)  The  variation 
of  n  -  1  with  wavelength  (dispersion)  has  been  studied  and  is  tabulated 
for  a  temperature  of  15°C,  which  must  be  a  convenient  laboratory  tempera¬ 
ture  in  Sweden.  The  following  data  applicable  at  15°C  are  obtained  from 
the  work  of  Edlen2: 

X( Angstroms )  5000  6000  7000  8000  9000  10000 

(n  -  1)  x  104  2.78  2.76  2.75  2.75  2.74  2.74 


At  the  more  comfortable  laboratory  temperature  of  20°C  the  following 


2.  B.  Edlen,  The  dispersion  of  standard  air,  J.  Opt.  Soc.  Amer.  43, 
339,  1953. 
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data  are  applicable  for  760  torr  pressure: 

Density:  p?Q  =  1.205  g  cnr3, 

Sound  velocity:  cJ0  -  3.44  x  104  cm  sec-1, 

Loschmid's  number  L  =  2.51  x  10 19  cm-3. 

The  value  of  n  -  1  varies  from  2.73  x  10"4  at  5000  A  to  2.69  x  10~4  at 
10000  A. 


The  attenuation  of  light  caused  by  scattering  on  the  molecules  of  air 
is  calculable  from  Rayleigh's  formula3 


o 


32a3(n-l)2 

3LX4 


(2.1) 


At  standard  temperature  and  pressure  we  obtain  the  following  values 
of  o: 

A(AngStrbms)  5000  6000  7000  8000  *000  10000 

a  x  107(cm’1)  1.70  0.81  0.44  0.25  0.16  0.10 

Both  (n-1)  and  L  are  proportional  to  p,  therefore  o  is  also  propor¬ 
tional  to  p.  At  constant  pressure  o  is  inversely  proportional  to  the 
absolute  temperature.  Therefore  the  values  of  o  at  t  =  20°C  and  760 
torr  are  about  7  percent  smaller  than  the  above  tabulated  values. 

3.  Absorption  Experiment  With  a  Steady  Laser  Beam. 

We  shall  analyze  the  problems  involved  in  the  measurement  of  «  by 
rowans  of  the  following  experiment:  A  certain  volume  cf  air  is  confined 

3.  J.W.  Strutt,  Lord  Rayleigh,  On  the  transmission  of  light  through  an 
atmosphere,  Phil.  Mag.  47_,  374,  1899;  Scient.  Papers  4,  397. 
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at  a  known  initial  pressure.  A  laser  beam  passes  through  the  confined 
volume  and  heats  the  enclosed  air  at  a  rate  proportional  to  a.  In  a 
steady  state,  the  heat  losses  of  the  system  balance  the  heat  gained  from 
the  laser.  The  increase  of  the  static  pressure  within  the  vessel  is  meas¬ 
ured  and  is  used  for  computing  the  rate  at  which  heat  is  added  tc  the  air. 

The  simplest  experimental  embodiment  of  this  idea  consists  of  a  cylin¬ 
drical  tube  kept  at  a  constant  temperature.  The  laser  beam  passes  along 
the  axis  of  the  cylinder  heating  a  smaller  cylindrical  volume  of  air  con¬ 
centric  with  the  confining  tube.  A  steady  state  is  established  when  the 
rate  at  which  heat  is  conducted  to  the  wails  is  equal  to  the  rate  at  which 
heat  is  absorbed  from  the  laser  beam. 

In  order  that  a  technique  of  this  type  yield  a  valid  measurement  of 
a,  certain  physical  conditions  must  be  satisfied  which  are.  .ct  easy  to 
achieve  experimentally.  First,  the  confining  tube  must  be  long  enough  so 
that  special  thermal  effects  a.t  the  ends  of  the  tube  do  not  significantly 
affect  t'.ie  pressure  of  the  confined  air.  Laser  light  is  brought  in  through 
solid  windows  whose  density  is  about  1000  times  that  of  the  confined  gas. 

It  is  net  unreasonable  to  assume  that  absorption  per  unit  volume  in  the 
windows  will  be  about  1000  times  as  great  as  absorption  in  the  gas. 

Second,  for  X  <  1.2u  more  laser  light  will  be  scattered  than  directly  ab¬ 
sorbed.  The  energy  of  the  scattered  radiation  must  be  removed,  it  must 
not  be  allowed  to  add  to  the  thermal  energy  of  the  confined  gas.  This 
removal  could  be  accomplished,  for  example,  by  making  the  inner  surface 
of  the  tube  black  and  by  making  the  tube  of  a  material  with  good  thermal 
conductivity. 
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Assuming  that  these  experimental  difficulties  are  mastered,  we  calcu¬ 
late  the  pressure  increase  dp  attained  in  the  vessel  in  the  steady  state 
as  a  result  of  energy  being  absorbed  from  the  laser  beam  at  a  constant 
rate  and  conducted  away  through  the  air  to  the  walls  of  the  vessel  which 

are  kept  at  a  constant  temperature  T  .  We  assume  a  uniform  circular  laser 

o 

beam  of  radius  a  in  a  vessel  of  radius  b.  Let  the  power  incident  in  the 
laser  beam  be  W,  then  the  intensity  of  the  beam  is  I  =  W/a2*.  Keat  is 
deposited  in  the  unit  volume  of  the  gas  at  the  rate  oJI,  where  J  =  0.239 
cal/ joule . 


I 

1 


i 
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Once  the  steady  state  is  reached,  no  mechanical  work  is  done.  There¬ 
fore  the  heat  balance  equation  is  obtained  by  equating  the  rate  at  which 
energy  is  absorbed  to  the  rate  at  which  it  is  conducted  away.  (It  is 
assumed  that  convection  is  neglected.)  The  rate  of  heat  generation  per 
unit  volume  is 

q'(r)  =  aJI  when  r  <  a, 

8  0  when  r  >  a. 

The  heat  balance  equation  for  a  volume  V  is  then 

/q'dv  =  -  K  ^VT-ndS,  (3.1) 

V  S 

where  dS  is  the  surface  element,  n  is  the  unit  normal  to  the  surface,  and 
K  the  thermal  conductivity.  It  follows  then  in  a  well-known  manner  that 
V2T  =  -  q'/K.  (3.2) 

Since  the  only  relevant  coordinate  in  this  problem  is  r,  the  distance 
from  the  cylinder  axis,  we  have 

Id  dT  q'(r) 

r  dr  r  dr  =  '  K  C3.3) 
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I 
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I 

I 

| 
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:  -1 

i  \ 

1  I 


This  differential  elation  is  to  be  solved  wich  the  boundary  condition 
T  =  Tc  for  r  =  b  and  the  solution  is  required  to  remain  finite  for  r  <  b. 
Although  che  solution  is  easily  found,  we  do  not  need  its  explicit  form 
because  we  are  on.y  interested  in  the  pressure  rise  in  the  vessel  and  this 
is  proportional  to  the  increase  of  the  average  temperature. 

In  place  of  T  we  introduce  the  dimensionless  variable 

V(r)  =  Sj T(T  -  To>-  (3.4) 

Equation  (3.3)  is  then  replaced  by 
Id  _  dv  .  /  . 

Fd?rdF- "  8(r)-  <3-s> 

where  g(r)  =  1  for  r  <  a  and  0  for  r  >  a.  The  boundary  condition  is 
v(b)  =  0.  Then  we  obtain  by  integration  of  (3.5) 

,  r 

r  dF  =  '  /  S^'Jr’dr'  =  -  f(r),  (3.6) 

o 

where  f(r)  =  j  r2  for  r  <  a,  and  f(r)  =  i  a2  for  r  >  a. 


The  average  increase  in  gas  temperature  over  the  entire  vessel  is 
proportional  to 


_  -  u  *  n  u  D 

v  =  7  /  /  /  v(r)rdrd6dz  =  K  J  v(r)rdr. 


b  2n  L 

Iff 

O  O  O  u  0 

Integrating  by  parts  and  noting  that  v(b)  =  0,  we  get 


(3.7) 


b  b  b 

/  v(r)rdr  =  -I  /  r2  p.  dr  =  A  /  f(r)rdr. 
o  2  q  dr  2  q 

b  a  b 

Since  /  f (r)rdr  =  /  i  r3  dr  +  /  *  a  rdr  =  [2  -  «i],  we  have 

o  o  2  a  2  8  b7 

v  =  —  ci  -  sLj. 

4  2b2 


(3.8) 
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Therefore  the  average  increase  cf  temperature  in  the  vessel  is 


dT  =  -J-I-d2  [i  . 

4K  2b2 


(3.9) 


The  final  expression  for  the  increase  of  pressure  in  the  chamber  is  ob¬ 
tained  from  (3.9)  by  using  the  gas  law  and  by  reintroducing  the  incident 
power  W  =  a2i»I. 


dp 


=  1°  (i 

"  T0  4*K 


(3.10) 


Ordinarily  the  expression  in  the  parentheses  is  nearly  equal  to  one. 
Assuming  a  laser  with  100  watts  useful  output,  an  initial  temperature 
Tc  =  300°K  and  a  =  10-9,  the  expression  in  front  of  the  parentheses  is 
1.1  x  10~7  atmospheres.  A  pressure  change  of  this  order  of  magnitude 
should  be  detectable  provided  one  can  keep  the  temperature  of  the  walls 
constant  to  within  0.3  x  lO-*4  °C.  It  might  be  quite  difficult  to 
accomplish  this  for  a  tube  as  long  as  1  meter. 


Dissipation  of  Heat  in  Air  by  Conduction. 

In  preparation  for  the  discussion  of  experiments  involving  the  heat¬ 
ing  of  air  colut  '.s  by  means  of  short  laser  pulses  we  shall  calculate  the 
rate  at  which  thermal  energy  initially  communicated  to  an  air  column  is 
dissipated  by  thermal  conduction  alone. 

The  calculation  is  simplified  by  neglecting  for  the  purpose  of  this 
calculation  the  aerodynamic  properties  of  air.  Thus  we  shall  be  concerned 
with  a  solid  substance  characterized  by  a  heat  conductivity  K,  density  p. 


i  i 
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and  specific  heat  cv. 


It  is  shown  in  textbooks4  that  the  temperature  v(x,y,z,t)  satisfies 
the  differential  equation 

c.i> 

provided  that  there  are  no  sources  or  sinks  of  heat  in  the  region  under 
consideration.  Here  k  =  K/pcv  is  the  diffusivity  of  the  substance. 


We  will  assume  a  geometry  of  cylindrical  symmetry  and  calculate  the 
cooling  of  an  infinite  cylinder  of  radius  a,  whose  temperature  at  the  time 
t  =  0  is  one  degree  above  the  temperature  of  its  surroundings.  All  space 
is  assumed  to  consist  of  the  same  material  characterized  by  the  constant 
k.  The  problem  then  requires  the  use  of  only  one  spatial  variable,  r, 
the  cylindrical  polar  coordinate.  The  differential  equation  becomes 


3v  . 
3t  * 


ill), 

r  3r  / 


and 


the  initial  condition 
f(r)  =  1 
=  0 


is  v(r,0)  =  f(r), 
for  r  <  a, 
for  r  >  a. 


where 


(4.2) 


As  long  as  no  walls  are  present,  no  boundary  conditions  are  required. 
This  differential  equation  is  readily  solved  by  the  Hankel  transformation 
technique^  and  the  following  result  is  obtained: 

4.  H.S.  Carslaw  and  J.C.  Jaeger,  Conduction  of  Heat  in  Solids,  Oxford,  1959 

2nd  ed.  p.  9.  ’ 

5.  Ref.  4,  p.  460. 
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v(r,t)  =  /  r0(o)oJ0(or)e“l<a2tdo, 


(4.3) 


where 


F0(a)  =  /  rJQ(or)f(r)dr. 
o 

In  the  present  case 

a 

Fq(°)  =  /  rJQ(or)dr  =  J^(oa). 


(4.4) 


(4.5) 


Hence  the  expression  for  the  temperature  distribution  is 


v(r,t)  =  a  /  J1(oa)J0(cr)e~<o2't  do. 


(4.6) 


It  is  interesting  to  calculate  the  variation  of  temperature  at  the  axis. 
From  (4.6)  we  have  for  r  =  0 


v(0,t)  =  a  /  J0(oa)e"ie°2t  do. 
o 

This  expression  is  changed  by  integration  by  parts  into 


v(0,t)  =  1  -  2xt  /  J0(oa)<ia-,ccj2t  do. 
o 

From  the  standard  formula6 


/  J0(oa)e_P2c?2  odo  =  — i—  e 
o  2p2 


(4.7) 


(4.8) 


(4.9) 


we  get 


v(0,t)  =  1  -  e 


'4<t 


(4.10) 


6.  G.N.  Watson,  Theory  of  Bessel  Functions,  Cambridge,  University  Press, 
1958,  p.  393. 
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Thus  the  temperature  at  the  center  line  decreases  with  t  according  to  the 
formula  1  -  e^/t,  where  T  =  a2/4<i  The  quantity  T  has  th#  dimension  of 

time.  It  is  the  time  during  which  the  temperature  decreases  by  the 
amount  1/e.  Although  we  are  not  dealing  here  with  a  simple  exponential 
decay  of  temperature,  it  is  still  reasonable  to  regard  t  as  a  thermal  re¬ 
laxation  time.  It  is  proportional  to  the  cross  sectional  area  of  the 
beam. 


Numerical  computation  of  the  thermal  relaxation  time  requires  the 
knowledge  of  *  for  air  in  the  applicable  temperature  range.  This  quan¬ 
tity  entered  first  in  equation  (4.1)  which  is  taken  ^rom  the  theory  of 
heat  conduction  applicable  to  solids.  When  this  theory  is  applied  to 

a  gas,  the  question  arises  what  kind  of  specific  heat  is  to  be  used  in 
the  equation 

x  -  K 

pF  (4.11) 

that  serves  as  a  definition  of  x2.  Since  we  are  contemplating  a  steady 
state  situation  in  which  no  mechanical  motion  takes  place,  no  work  is 
being  done  and  the  entire  heat  energy  that  enters  a  volume  of  gas  is 
converted  into  internal  energy.  The  specific  heat  at  constant  volume 
is  the  one  which  is  applicable  to  such  a  situation.  For  air  at  standard 
pressure  and  temperature,  x  =  0.255  cm2  see"*.  Although  x  depends  on  the 
t«  .iperature,  this  value  of  x  is  probably  a  fair  approximation  to  the 
correct  value  applicable  at  20°C.  With  the  value  x  =  0.255  cm2  sec-1 
we  obtain  the  following  approximate  values  of  relaxation  times  and  as.oci- 
ated  frequencies  applicable  to  beams  of  convenient  cross  sections. 
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Beam  diam.  (2a)  mm 

0.4 

0.6 

0.6 

1.0 

1.2 

1.4 

1.6 

1.8 

2.0 

Relax,  time  (t)  msec 

0.392 

0.882 

1.57 

2.45 

3.53 

4.80 

6.27 

7.94 

9.80 

Frequency  (t_1)  Hz 

2550 

1134 

637 

408 

283 

208 

160 

126 

102 

In  a  dynamic  situation, 

when 

heat  is 

added 

to  a 

column 

of  gas,  not 

only  heat  conduction  takes  place  but  the  gas  is  set  in  motion.  This  trans¬ 
ient  phenomenon  was  not  taken  into  account  so  far.  In  what  follows  later, 
the  dynamic  phenomenon  is  of  great  importance  and  the  results  of  this  section 
are  used  for  guidance  concerning  the  situations  when  the  effect  of  heat  con¬ 
duction  may  be  neglected.  Clearly  one  can  only  justify  neglecting  heat  con¬ 
duction  if  the  dynamic  phenomenon  (oscillation)  occurs  so  fast  that  only 
negligible  cooling  takes  place  during  one  cycle.  Such  will  be  the  case 
when  the  frequency  of  the  sound  wave  generated  exceeds  the  appropriate 
characteristic  frequency  t'  1  =  4</a2. 

When  the  air  column  is  not  free,  as  heretofore  assumed,  but  confined  in 
a  pipe  kept  at  a  constant  temperature,  the  method  of  calculation  is  different. 
The  differential  equation  (4.1)  is  applicable  in  this  case,  but  a  boundary 
condition,  absent  in  the  case  of  the  free  air  column,  must  be  added.  For  a 
cylindrical  vessel  of  inner  radius  b,  the  boundary  condition  is  v(b,t)  =  0, 
provided  that  the  vessel  is  kept  at  a  constant  temperature  T0.  (Here  v 
denotes  the  excess  over  this  temperature.)  The  solution  of  the  heat  con¬ 
duction  problem  for  a  domain  with  a  boundary  is  well  known7.  It  is  obtained 
by  expanding  the  function  to  be  determined  in  terms  of  the  eigenfunctions 
of  the  differential  equation  (4,1)  with  the  given  boundary  conditions. 

In  the  case  of  heat  conduction  within  a  cylinder  of  radius  b  the  result  is 

Ref.  4,  pp.  198-199. 
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(4.12) 


v(r,t)  =  [  AnJ0(anr)e-<<sn:, 
n=l 

where  the  numbers  an  are  determined  from  the  boundary  condition  JQ ( b )  =  0, 
and  the  coefficients  Af)  are  calculated  from  the  expansion 

m 

s  l  AnJo(<*nr),  (4.13) 

n=l 

which  is  the  explicit  form  of  the  initial  condition.  The  An»s  are  deter¬ 
mined  by  the  routine  procedure  applicable  to  expansions  according  to  or¬ 
thogonal  functions.  Their  values  depend  on  a  and  b.  In  the  general  case 
the  value  of  AR  does  not  drop  off  rapidly  with  increasing  n.  We  must 
therefore  take  into  account  a  large  number  of  terms  in  (4.12).  Each  term 
decays  in  time  according  to  a  different  exponential  rate  kci2.  The  time 
constants  associated  with  these  modes  are  Ua2)-1,  and  an  is  deter¬ 
mined  from  the  equation 

“nb  s  xn»  (4.14) 


where  x  is  the  n-th 
n 

root  of  the  equation  JQ(x)  =  0. 

The  values  of  the 

first  five  x's  are: 

2.405,  5.520,  8.654,  11.792  and 

14. 913. 8  The  time 

constants  are 

b2 

Tn 

*Xn  ' 

(4.15) 

We  note  that  the  role  of  the  Tn's  is  not  the  same  as  that  of  t  in 
the  case  of  the  free  cylinder  because  in  the  present  case  the  time  de¬ 
pendence  of  a  mode  is  of  the  form  exp(-t/tn)  and  not  1  -  exp(-T/t). 

By  means  of  a  physical  argument  we  can  readily  infer  that  the  confined 

8.  A.  Gray  and  G.B.  Mathews,  Bessel  Functions,  2nd  ed.,  Dover  Publ.,  New 
York  1966,  p.  300. 
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cylinder  will  cool  more  rapidly  than  the  free  cylinder  and  that  the  re¬ 
laxation  time  calculated  for  the  free  cylinder  will  be  an  approximation 
of  the  effective  relaxation  time  of  a  confined  cylinder  for  b  >>  a. 

5,  Pressure  Rise  in  a  Tube  Caused  by  a  Single  Laser  Pulse. 

One  might  attempt  to  confine  air  in  a  cylindrical  tube,  shoot  a  laser 
beam  through  and  measure  the  rise  in  pressure  in  the  tube.  The  expected 
increase  in  pressure  is  proportional  to  the  energy  deposited  in  the  air 
by  the  laser  beam.  Hence  one  expects  thv.t  the  absorption  coefficient  a 
might  be  measured  in  such  an  experiment. 

It  must  be  remembered  that  the  walls  of  the  tube  holding  the  air  are 
solid  and  have  a  much  larger  heat  capacitance  than  the  confined  air.  So 
eventually  as  thermal  equilibrium  is  approached,  the  thermal  energy  de¬ 
posited  in  the  air  is  conducted  to  the  walls,  temperature  and  pressure  in 
the  vessel  return  to  their  values  at  the  beginning  of  the  experiment.  The 
approach  to  thermal  equilibrium  may  be  thought  of  as  consisting  of  a  fast 
process  of  pressure  equalization  followed  by  a  slower  process  of  tempera¬ 
ture  equalization  by  heat  conduction.  The  pressure  equalization  may  be 
characterized  by  the  time  constant  2b/c  where  b  is  the  radius  of  the  tube 
and  c  the  velocity  of  the  sound.  The  rate  at  which  heat  conduction  leads 
to  thermal  equilibrium  must  be  calculated  in  the  manner  described  in  Section 
4.  The  time  constant  t  =  b2/4<  certainly  reveals  the  correct  order  of  mag¬ 
nitude.  The  measurement  of  pressure  increase  should  take  place  after  the 
elapse  of  the  time  2b/c  following  the  passage  of  the  laser  beam  but  consid¬ 
erably  before  the  elapse  of  the  time  t.  For  a  tube  of  1  cm  radius, 

2b/c  =  6  x  10"5  sec  and  t  =  1  sec. 


16  - 


The  laser  beam  will  generally  be  much  narrower  than  the  tube.  Let 
the  radius  of  the  beam  cross  section  be  i,  and  let  us  assume  that  we  are 
dealing  with  a  very  short  laser  pulse.  Then  at  the  end  of  the  pulse  we 
have  an  air  cylinder  of  radius  a  within  which  the  pressure  has  been  raised 
to  P  ♦  dp  and  the  temperature  to  T  +  dT,  while  outside  this  cylinder  the 
pressure  is  p  and  the  temperature  is  T.  According  to  the  standard  equation 
of  heat  balance 

cvpoa2ffdT  =  aJQ*  (5.1) 

where  Q  is  the  energy  of  the  laser  pulse  in  joules,  J  is  the  conversion 
factor  0.239  cal/ joule,  and  pc  is  the  density  of  the  air.  From  the  ideal 
gas  law  and  equation  (5.1),  we  get 


dp  _  ctjQ 
p  cvpoa2irT  ' 


(5.2) 


On  substituting  the  appropriate  physical  constants  of  air,  we  get 
for  a  =  0.1  cm  and  T  =  300°K: 

^  =  1.14  x  102  oQ.  (5.3) 

Thus  a  laser  pulse  of  1  joule  and  an  a  =  10~9  cm~^  would  give  a  relative 
pressure  increase  of  one  part  in  ten  million.  It  seems  then  that,  using 
a  powerful  laser  of  several  joules  output,  an  a  of  the  order  of  10~9  cm-1 
would  he  just  barely  detectable.  This  would  be  the  case  if  the  entire  con¬ 
fined  space  would  be  within  the  beam.  Actually  the  beam  diameter  must  be 
considerably  smaller  than  the  diameter  of  the  tube  that  keeps  .he  air  con¬ 
fined.  Consequently  the  pressure  rise  observable  in  the  tube  is  much  less 
than  the  pressure  rise  in  the  cylinder  swept  out  by  the  beam. 


An  estimate  of  the  pressure  rise  anticipated  in  the  tube  may  be  ob¬ 
tained  by  ascribing  the  entire  gas  in  the  tube  some  soi't  of  average  tem¬ 
perature  T.  The  heat  balance  equation  (5.1)  is  then  replaced  by 

cvp0b2ndf  =  aJQ.  (5.4) 

This  equation  differs  from  (5.i)  insofar  as  b2  takes  now  the  place  of  a2. 
Consequently  the  pressure  rise  calculated  from  (5.4)  will  be  reduced  in 
proportion  to  the  ratio  of  the  cross  sectional  ratio  of  the  beam  to  the 
cross  sectional  area  of  the  tube.  If,  for  example,  a  =  0.1  cm  and  b  =  1.0  cm 
ti  e  value  of  cip/p  is  reduced  by  a  factor  of  100,  thus  driving  it  completely 
beyond  reach  for  measurement . 


It  could  be  objected  that  the  above  estimate  of  dp/p  is  based  on  a 
dubious  assumption  of  an  average  temperature  and  the  application  of  thermo¬ 
dynamics  to  non-equilibrium  states.  The  deduction  may  be  replaced  by  the 
following  one:  After  irradiation,  the  pressure  of  the  air  within  the  beam 
is  p^,  its  volume  is  v^.  In  the  rest  of  the  vessel  the  pressure  is  pQ. 

The  volume  of  the  rest  is  V2*  Both  volumes,  change  during  the  process  of 
pressure  eque lizatics;  which  results  in  a  final  pressure  P.  Assuming  no 
heat  exchange, 

vx  =  pl/Y  v|  and  v2  =  P*^  v2* 


Moreover,  the  total  volume  does  not  change,  therefore  vi  +  v2  =  vl  +  v2* 
Hence 


,1/Y  -  pl 


1/Y 


'1  ^  pJ/Y 


(5.5) 


V1  ♦  v2 


and 


,1/Y  _  1/Y 


V1  +  v2 
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(5.6) 


From  this  last  equation  it  also  follows  that  in  first  approximation 


6.  Generation  of  Sound  by  a  Laser  Pulse. 

A  laser  pulse  is  sent  through  a  gas  which  has  a  small  absorption  co¬ 
efficient  for  the  radiation  sent  through  it.  As  a  corsequence  of  absorption 
a  column  of  gas  is  warmed  and  the  surrounding  gas  is  set  in  motion.  It  is 
desired  to  determine  the  sound  energy  generated  by  the  pulse  and  to  determine 
the  spectral  distribution  of  the  sound. 

The  solution  of  the  problem  is  expected  to  contain  the  physical  para¬ 
meters  of  the  gas  which  normally  determine  the  velocity  of  propagation  of 
sound  in  the  gas,  the  rate  at  which  radiated  energy  is  absorbed  in  the  gas, 
the  length  of  the  laser  pulse  and  the  cross-section  of  the  laser  beam. 

It  will  be  assumed  that  the  gas  is  an  ideal  gas,  that  the  laser  beam 
is  of  circular  cross-section,  that  the  rate  of  energy  dissipation  and  the 
pulse  length  are  such  that  the  resulting  temperature,  pressure  and  density 
changes  may  be  treated  as  differentials  of  the  first  order  so  that  the 
usual  approximations  employed  in  the  theory  of  sound  retain  their  validity. 

The  effects  of  heat  conduction  will  be  neglected;  a  simplification 
which  may  only  be  questionable  when  it  is  applied  to  low  frequencies.  The 
low  frequency  limit  of  this  approximation  is  the  characteristic  frequency 
associated  with  the  cooling  of  the  heated  column  by  conduction.  This  was 
calculated  in  Section  4, 


The  propagation  of  sound 


waves  is  described  by  a  solution  of  the 


wave  equation 


^f-c2'2*  =  3  (6.1) 

which  is  satisfied  by  the  velocity  potential  41  as  well  as  by  the  pressure 
p  and  the  condensation  s  in  the  gas.  The  differential  equation  (5.1)  is  a 
homogeneous  one,  the  solution  $(x,y,z;t)  is  determined  if  one  knows  the 
value  of  +  and  3$/3t  at  the  time  t  =  0. 


Th.  physical  problem  to  be  solved  here  must  be  governed  by  the  dif¬ 
ferential  equation  (6.1)  in  the  region  where  no  heating  takes  place.  In 
the  region  traversed  by  the  laser  beam  we  expect  to  obtain  an  inhomogeneous 
differential  equation  whose  solution  is  to  be  obtained  for  a  given  right 
side.  Moreover  the  solution  is  to  vanish  for  t  =  0;  it  shou i i  represent 
an  outgoing  wave  away  from  the  region  traversed  by  the  laser  beam  and  is 
to  satisfy  the  usual  continuity  conditions. 


The  — -rting  point  in  acoustics  is  the  continuity  equation  for  the 

fluid 


7-pV 


0 


(6.2) 


and  the  equation  of  motion  for  the  fluid  particles: 
dV 

p  dt  =  ‘  7p'  (6.3) 

Here  the  vector  V  denotes  the  velocity  of  the  fluid,  p  its  pressure,  p 
its  density.  All  these  quantities  are  functions  of  position  and  time. 
It  is  customary  to  restrict  the  problem  to  small  velocities  and  small 
displacements  from  static  or  average  pressure  and  density  pQ  and  pQ  . 
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With  these  restrictions  equations  (5.2)  and  (6.3)  are  replaced  by 


and 


<jp  _ 

at  +  Po7*v  =  0 


3V 

Po  ^  t  7p  =  °. 


We  introduce  the  velocity  potential  $  defined  by 
V  =  -V*, 


then  (6.5)  becomes 


or 


•PoI7*  ♦  VP  =  O' 


7(|i  -  £_)  =  0. 

at  pQ 


(6.4) 


(6.5) 


(6.6 ' 


(6.7) 


(6.8) 


We  write  p  =  p  +  p\  where  p  is  a  constant  and  p'  is  a  variable 


quantity 


small  compared  to  pQ.  Then 


3£  _  El 

3-  "  Po  ’ 


(6.9) 


This  last  equation  relates  the  pressure  variation  to  the  velocity  potential. 


Now  the  continuity  equation  (6.4)  is  rewritten  by  introducing  the  con¬ 
densation  s  defined  as  follows 


or 


P  =  pQ(l  +  s), 


P  -  Pt 


Po 


(6.10) 


(6.11) 


Then  (6.4)  becomes 

3s 

Jt 

or 


+  v-v  =  o. 


(6.12) 
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(6.13) 


3s  . 

at 


?2<j>. 


In  the  classical  case  the  equation  of  sound  propagation  |_Eq.  (6.1)] 
is  found  by  assuming  that  the  pressure  variation  p'  is  p*  portional  to 
the  condensation  s.  The  factor  of  proportionality  is  found  Hy  thermo¬ 
dynamic  argument  to  be  equal  to  ypD.  With  this  result  one  gets 


P*  =  YP0S- 


(6.14) 


Hence  from  (6.9) 

32<f>  yp0  3s 
3t2  pQ  3t 


Finally 


32$ 

3^ 


(6.15) 


(6.16) 


which  is  id  -tical  with  equation  (6.1).  Its  derivation  is  based  on  the 
assumption  that  no  thermal  exchange  takes  place. 


When  the  gas  is  locally  heated  equation  (6.14)  is  no  longer  valid. 
We  must  calculate  the  change  in  pressure  when  both  the  volume  is  changed 
and  . .  certain  amount  of  heat  is  added.  If  the  heat  is  added  reversibly, 
we  write  dQ  =  TdS.  Then 

dp  =  (|E)  dV  +  (||)  dS.  (6.17) 

dV  g  do  y 

Equation  (6.14)  is  a  special  case  of  (6.17)  for  dS  =  0. 


In  the  general  case  the  coefficients  of  the  differentials  in  (6.17) 
are  evaluated  as  follows:  The  entropy  change  of  1  mole  of  an  ideal  gas 
may  be  written  in  the  following  two  forms: 9 

9.  M.W.  Zemansky:  Heat  and  Thermodynamics,  McGraw-Hill,  New  York,  .1943, 
Sections  13.2  and  13.3. 
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ds  =  cpf-^- 


(6.18) 


=  dV. 

(6.19) 

Upon  elimination  of  dT  from  (6.18)  and  (6.19)  we  obtain 

-  dp  =  (y  -  l)dS  -  y-  dV. 

P  T 

(6.20) 

Hence  from  (6.20)  and  the  gas  law  pV  =  RT  it  follows  that 

dp  =  -  yjl  dV  +  -Y-  ~  1iTdS, 

(6.2i) 

and  therefore 

<i£)  =  .  m 

3V'S  V  ’ 

(6.22) 

(ie.)  s  y  -  i  T 

l3S  v  V  T* 

(6.23) 

We  wish  to  express  the  pressure  change  p'  in  terms  of  the  condensation 

s  and  the  heat  added  per  unit  volume  q.  These  quantities  are  related  to 
dV  and  dS  as  follows 


Thus 


P'  =  *PoS  +  <Y  -  Dq.  (6.24) 

Equation  (6.24)  must  be  used  in  place  of  (6.14)  to  obtain  the  equations 
valid  in  the  case  of  the  heated  gas  column.  Differentiating  (6.9)  we 
have 

3p^  YPo  ?s  Y  -  1  3q 

3t2  "  po  3t  =p7”3t  +  _pl  it*  (6,25> 

On  combining  (6.25)  with  the  continuity  equation  (6.13),  we  get 
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(6.26) 


32$  YP0  Y  -  1  dq 

77  s  —  7  *  + - • 

3t2  pQ  pg  at 

In  the  absence  of  heat  transfer  by  conduction  3q/9t  is  equal  to  q' ,  the 
rate  of  heat  generation  per  unit  volume.  When  thermal  conduction  is  taken 
into  account  then 


+  KV2T, 


(6.27) 


in  accordance  with  the  heat  balance  equations  (3.1)  and  (3.2).  Thus, 
when  conduction  is  neglected,  we  have 

.  C2V2$  =  ,  (6.28) 

3t2  pQ 

where  «  =  /ypQ/po  is  the  velocity  of  sound.  It  is  apparent  from  this 
equation  that  of  the  energy  communicated  to  the  gas  only  the  fraction 
y  -  1  is  converted  into  an  acoustic  (hydrodynamic)  wave.  In  the  case  of 
diatomic  gases,  this  fraction  is  0.40.  Given  a  pulsed  laser  of  such  total 
output  that  Q  is  the  energy  incident  upon  the  test  gas  d  ing  a  period  t, 
the  rate  of  heat  deposition  per  unit  volume  is 


q'  =  oO/At,  (6.29) 

where  A  is  the  cross  sectional  area  of  the  laser  beam  and  a  the  absorption 
coefficient.  The  differential  equation  (6.28)  may  then  put  in  the  form 


where 


V2$ 


~  ~  -  g(x,y,z,t), 

c‘  3t2 


(6.30) 


8 


Y  "  1  —  f(x,y,z,t). 
PQc2  At 


(6.31) 


Here  f  denotes  a  function  which  is  1  within  the  laser  beam  during  the 
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duration  of  the  pulse  and  0  at  every  other  point  in  space  or  outside  the 
time  interval  0  <  t  <  t. 

The  solution  of  this  differential  equation  depends  on  the  initial  and 
boundary  conditions.  We  assume  the  homogeneous  initial  conditions 
$(x,y,z,0)  8  0  and  $t(x,y,z,0)  =  0  which  correspond  to  no  disturbance  at 
t  =  0.  We  then  have  a  standard  radiation  problem  solvable  in  terms  of  the 
retarded  potentials.  In  the  absence  of  bounda' y  conditions;  i.e.  for  a 
homogeneous  medium  of  infinite  extent,  the  solution  of  (6.30)  for  homo¬ 
geneous  boundary  conditions  is  given  by  the  formula:10 

$(x,y,z,t)  =  ~  Iff  g(£,n,S,t-r/c)r-1  d£dndc,  (6.32) 

where  r  =  /(x-£)z  +  (y-n)2  ♦  (z-?)2.  The  function  g  is  a  constant  times 
a  step  function.  We  write  Eq.  (6.31)  in  the  form 

g  *  ~  f(x,y,z,t),  (6.33) 

T1 

where 

Ap0c2t 

T1  (y-DctQ  (6.34) 

is  a  constant  that  has  the  dimension  of  time.  Then 

»(x,y,z,t)  *  JL-  ///  *< *,.*[•  C**‘r/c)  dCdndc.  (6.35) 

The  integration  is  to  be  extended  over  the  volume  of  space  swept  by  the 
laser  beam.  Its  result  depends  only  on  the  beam  cross  section  and  the 
length  of  the  pulse. 

10.  A.N.  Tikhonov  and  A. A.  Samarskii,  Equations  of  Mathematical  Physics, 
MacMillan,  New  York,  1963,  pp.  461-465. 
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To  evaluate  the  integral  in  (6.35)  we  place  the  z-axis  along  the 
direction  of  the  propagation  o  the  laser  beam  and  we  introduce  cylin¬ 
drical  coordinates  R,  and  z.  Assuming  circular  symmetry,  we  obtain 
readily 


♦<R.t)  =  -A_  / 
4irTl  g 


d; 

*k2+(z-c)2  ’ 


(6.36) 


where  the  region  of  integration  is  confined  to  such  values  of  £  ior 
which  0  <  t  -  r/c  <  t ,  or  equivalently 

t  -  t  <  r/c  <  t.  (6.37) 

Equation  (6.36)  is  an  approximation  valid  when  the  distance  R  from  the 
center  of  the  laser  beam  is  much  larger  than  the  diameter  of  the  laser 
beam,  i.e.  when  one  may  replace  (x-£)2  +  (y-n)2  by  x2  +  y2  =  R2  for 
every  point  (£*n,c)  within  the  laser  beam.  We  introduce  the  new  variable 
C'  =  z-5,  then  equation  (6.36)  becomes 


*(R,t)  = 


A 

4irtjL 


/r2+c'  2 


(6.38) 


where  the  range  of  integration  over  is  adjusted  to  be  consistent  with 
(6.37).  This  consistency  moans  that 

c2(t-x)2  -  R2  <  c'2  $  c2t2  -  R2.  (6.39) 


The  integration  in  (6.38)  is  symmetric  about  c'  =0,  therefore  we 


write 


*(R,t)  =  -i-  / 
2nt,  ; 


U2 


dC 


(6.40) 


where  and 
now  distinguish 
a )  When  R 


are  non-negative  limits  consistent  with  (6.39).  We  can 
three  regions: 

>  ct,  the  disturbance  has  not  reached  the  obse^  c-r>. 


A 
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=  0  and  $  =  0. 

b)  When  ct  >  R  >  c(t-T),  the  head  of  the  disturbance  has  reached 
the  observer,  but  its  tail  has  not.  In  this  case  =  0,  ^  0 

and 


x  _  A  ,  ct+»^t 

4  -  2^1  1o* - r 


(6.41) 


c)  When  R  <  c(t-x),  the  tail  of  the  disturbance  has  begun  to  reach 
the  observer.  In  this  case  u  $  0  and  u  I  o  and 


(log  (ct+/c^t5-Rz)  -  log  [c(t-t)  +  /cJ(t-?)z-Rz]  (6.42) 


$(cR 1  ,t )  =  {log  (t+/tz-R,z)  -  log  [t-T+/(t-T)z-R,z]}  (6.44) 

for  R*  <  t  -  t.  The  graph  of  $  as  a  function  of  the  variable  t  _s  sketched 
in  Figure  1  for  two  different  values  of  R'.  The  spectral  distribution  of 
the  function  4  is  clearly  not  the  same  as  that  of  the  source.  The  function 
♦  is  not  a  square  wave.  Its  shape  depends  on  the  distance  R’.  It  is  ex¬ 
pected  intuitively  that  ar.  observer  near  the  z-axis  would  observe  something 
more  resembling  a  square  wave  than  an  observer  located  farther  away  in  a 
region  where  the  sharp  corners  of  the  excitation  function  are  washed  out. 
The  spectral  distribution  is  calculated  as  follows: 
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The  time  derivative  o  $  has  a  very  simple  form.  In  fact. 


4*  =  =  0  for  R'  >  t, 

<*t 


2ittl  /£7_r 


— —  „■  for  t  >  R'  >  t-t, 

-R  ' 1 


(6.45) 


A 

2itt 


1  [  */t^-R  2/Tt-r)2-R,2j 


for  t-t  >  R' . 


The  following  equation  is  well  known: 


/  Jc(bx)  sin  axdx  =•( 
o 


for  a  <  b, 


(6.46) 


for  a  >  b. 


Replacing  b  by  R'  and  a  by  t  we  have 

m 

/  Jq(R'x)  sin  txdx  *  C  (6.47) 

o 

for  0  *  t  <  R'.  We  shall  restrict  our  observations  to  the  region  R'  > 
which  is  the  region  of  the  space  outside  of  the  cylinder  reached  by  the 
leading  edge  at  the  disturbance  when  the  pulse  ends.  We  integrate  (6.47) 
from  t  :  0  to  t  =  t  <  R'  and  obtain 


-J0(R'x) 


(1  -  cos  tx)  dx  -  0. 


(6.48) 


From  (6.45)  and  (6.46)  we  find  that  may  be  written  as  follows 


/  J0(Rfx)  sin  txdx  -  /  J0(R'x)  sin  (t-t)xdxl.  (6.49) 

o  o  J 

Integrate  according  to  t  from  t  =  0  to  t  =  s  and  note  that  ♦(cR',0)  =  0. 


11  _  A 
3t  2»ta 
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Then 


$(cR;-5)  =  _A_  j  J  (R'x)  C°s(t-T)x  _  COS  tx 

Out  L  x 


(6.50) 


Although  J0(R'x)  £Eiii  is  not  integrable  from  0  to  »,  the  integral  in 
(G.50)  is  convergent  because  the  expression  in  the  brackets  is 


cos  (s-t)x  -  1  ^  1  -  COS  TX  ^  1  -  COS  sx 
}:  X  X 


and  the  product  of  each  of  these  functions  with  J0(R'x)  is  integrable. 
In  fact  the  middle  term  contributes  nothing  (see  eq.  6.48).  Then 

^cR',s)  =  27fT-  ^  ^ -  £Cos  (s_T>x  -  cos  sx]  dx.  (6.51) 


Hence,  using  an  elementary  trigonometric  identity,  we  get 


♦('R'  ,s)  =  /"£ocp0  sln  jx  sin 


Now  let  t  =  s  -  y  and  x  =  w,  then 

*(CR'  /*  sin  fi  sin  »t 


(6.53) 


We  have  represented  *(cR\tf|)  as  a  Fourier  integral.  The  Fourier  sine 


transform  of  <t>  is  thus 


♦(»)  =  —  Jn(R'u))  £ln-  UT/2 


*T1  ° 


(6.54) 


Here  we  define  the  sine  transform  by  the  relations 


F(w)  =  1  /  f(t)  sin  ut  dt,  and  f(t)  =  /  FU)  sin  cat  dw. 
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The  spectral  tensity  of  <J>(cR',t)  is  therefore 


,)2  =  j  2(r’u,)  \2  _ii_ . 

°  \  T.  I  *2  2 


(6.55) 


It  is  shown  by  elementary  calculations  that  the  spectral  density  function 
of  a  square  wave  'f  unit  amplitude  and  of  duration  t  is 


.u))2  = 

s 


(6.56) 


Comparisor  of  (6.55)  and  (6.56)  shows  that  in  our  problei  the  spectrum  of 

the  square  wave  gets  distorted  by  the  factor  J  (R'w).  The  distortion  is 

o 

small  rear  the  source  where  Jq  %  1.  For  R’oi  %  2.4,  however,  Jq  =  0.  This 
means  that  at  a  given  distance  R  the  frequency  for  which  R  =  2.4c  will  be 
missing.  Oxher  missing  frequencies  correspond  to  other  zeros  of  JQ.  The 
reduction  ir  the  spectral  intensity  due  to  the  J0(R'u))2  factor  can  be  cal¬ 
culated  from  the  following  table  which  gives  the  reduction  from  the  inten¬ 
sity  of  a  square  pulse.  The  entry  item  of  the  table  is  the  product  Rv  in 
cm  times  cycles  per  second.  Assumed  velocity  c  =  3.3  x  104  cm/sec. 

vR(cm/sec)  500  1000  1500  2000  2500  3000  3500  4000  4500 

j  2f A  o  qg6  o.980  0.960  0.929  0.891  0.846  0.794  0.741  0.680 

o  \  c  / 

It  is  clear  that  only  for  large  distances  or  high  frequencies  will 
the  Bessel  factor  affect  significantly  the  spectral  distribution  of  the 
rgy.  The  main  determining  factor  at  small  oistances  and  moderate 
frecicncies  is  the  spectral  distribv  ion  of  the  pulse.  This  distribution 
is  characterized  by  the  function 


iillilillilliriliillililiiilliiilillliriiiliiiilil^  . .  . . . .  — 


which  is  well  known  from  the  theory  of  diffraction.  It  is  sketched  in 
Fig.  2.  The  largest  part  of  the  energy  is  contained  in  the  spectral  region 
between  0  and  uit  =  2tj.  Thus,  in  terms  of  frequency,  the  interval  0  <  v  <  1/t 
contains  most  of  the  energy.  Moreover,  it  is  manifest  from  Fig.  2  that  the 
spectral  distribution  is  rather  flat  in  the  region  in  which  most  of  the 
energy  is  contained.  A  pulse  of  this  type  is  not  suitable  as  a  source  for 
a  signal  which  is  tc  be  detected  by  a  narrow-band  detector. 

The  question  arises,  v;h*t  happens  when  a  pulse  of  the  type  described  is 
repeated  periodically  at  equal  time  intervals  p.  In  this  case  the  Fourier 
integral  is  replaced  by  a  Fourier  series  and  the  continuous  spectral  distri¬ 
bution  shown  in  Fig.  2  is  replaced  by  a  discrete  spectrum  consisting  of 
equidistant  lines  with  the  separation  Av  =  1/p.  The  relation  of  this  dis¬ 
crete  spectrum  to  the  continuous  spectrum  of  the  non-periodic  (single  pulse) 
case  is  illustrated  in  Fig.  3.  Again  the  energy  is  generally  broadly  dis¬ 
tributed  and  only  for  p  %  2t  does  a  significant  spectral  concentration  of 
energy  take  place. 

Note:  It  is  possible  to  remove  the  restriction  that  the  observation 
should  be  confined  to  the  region  R'  >  t.  This  assumption  was  made  following 
equation  (6.47).  In  tne  general  case  we  cannot  assert  that  the  left  hand 
side  of  equation  (6.48)  is  zero,  but  we  can  assert  that  it  is  a  function  of 
R'  and  t  alone.  Let  it  be  denoted  by  <KR',t).  Then  we  can' carry  out  the 
calculations  that  lead  tc  equations  (6.52)  and  (6.53)  provided  we  replace 
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{(cR'.s)  by  $(cR',s)  +  ii»(R’,r),  where  the  function  v  is  independent  ol  t. 
Consequently  the  spectrur.  of  $(cR',t)  will  be  tne  same  as  before  since  the 

time  independent  term  can  not  affect  for  u  4  0. 

7.  The  Acoustic  Prcnlem  in  a  Confined  Region. 

In  the  preceeding  section  we  analyzed  the  generation  of  a  sound  pulse 
from  the  thermal  energy  ceposited  in  a  gas  by  laser  oulses.  No  boundary 
conditions  were  imposed,  thus  implying  that  the  gas  was  not  confined.  If 
the  gas  were  not  confined,  or  if  the  confining  walls  were  far  removed  from 
the  path  of  the  beam,  acoustic  energy  would  radiate  in  all  directions.  The 
intensity  of  the  acoustic  signal  would  then  decrease  as  the  reciprocal  of 
the  distance  from  the  center  line  of  the  generating  beam.  The  detection  or 
measurement  of  such  a  signal  depends  to  a  large  extent  on  the  experimental 

arrangement  used  to  concentrate  a  significant  fraction  of  the  signal  gener¬ 

ated  on  a  suitable  detector  capable  of  converting  the  incident  acoustic 
signal  into  an  electric  signal.  It  would  require  a  rather  sophisticated, 
well-designed  system  to  utilize  as  much  as  20%  of  the  generated  acoustic 
energy.  Even  if  much  of  the  energy  radiated  out  is  concentrated  on  the 
detector,  by  means  of  reflectors,  for  example,  a  serious  degiadation  of 
the  detector  system  will  result  from  the  fact  that  the  generated  signal 
has  a  rather  broad  spectral  distribution. 

Now  we  shall  examine  the  generation  of  sound  in  a  confined  system,  a 
cylinder  concentric  with  the  laser  beam.  It  is  expected  that  in  such  a 
system  the  efficiency  of  the  utilization  of  the  sound  generated  '’ill  be 
higher,  and  one  migh*  hope  to  make  use  of  the  resonances  of  the  confining 
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tube  to  enhance  detection.  The  following  questions  are  worth  examining: 

1.  How  does  the  confining  tube  affect  the  generation  and  propagation 
of  the  sound  wave? 

2.  Is  it  possible  to  enhance  the  measurement  of  the  sound  wave  gener¬ 
ated  by  proper  choice  of  the  geometry  of  the  confining  vessel? 

The  starting  point  of  the  analysis  is  the  differential  equation  (6.30) 

?24>  "  =  -  g(x,y,z,t)  (7.1) 

c2  3t2 

with  the  initial  conditions 

♦(x,y,z,0)  =  0;  $t(x,y,z,0)  =  0.  (7.2) 

The  meaning  of  the  variables  was  explained  in  section  6,  where  it  was  ilso 
noted  that  g  is  zero  except  in  the  part  of  the  space  swept  out  by  the  laser 
beam.  In  this  cylindrical  region,  g  has  a  constant  value  during  the  time 
interval  0  <  t  .£  t;  it  is  0  when  t  >  t.  We  shall  simplify  the  computations 
by  normalizing  g  so  that  is  has  the  values  0  and  1.  Since  g  has  the  dimen¬ 
sion  of  reciprocal  time,  it  is  more  correct  to  say  that  we  set  the  peak 
value  of  g  at  1  sec-1. 

The  presence  of  the  walls  of  the  tube  requires  the  introduction  of  a 
boundary  condition  which  expresses  the  fact  the  normal  component  of  the 
fluid  velocity  at  the  walls  is  zero.  This  means 

n*V$  =  C.  (7.3) 

at  all  boundaries.  We  shall  seek  a  solution  of  (7.1)  in  the  region 
x2  +  y2  s  b2  for  t  i  0.  A  useful  technique  for  the  solution  of  such  a 
problem  is  the  eigenfunction  expansion.  This  technique  relies  on  the  deter¬ 
mination  first  of  the  simple  harmonic  solutions  of  the  homogeneous  equation 
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obtained  from  (7.1)  with  g  i  0.  It  is  well  known  that,  if  we  set 

$  =  e_iwt  i|i(x,y,z),  (7.4) 

then  4 1  satisfies  the  homogeneous  Helmholtz  equation 

V2*  +  k2*  =  0  (7.5) 

which  is  obtained  from  (7.1)  by  setting  g  =  0,  substituting  (7.4)  and 
k  =  w/c.  It  is  also  known  that,  given  a  finite  region  9?  and  g;ven 
homogeneous  boundary  conditions,  non-vanishing  solutions  of  (7.4)  are 
possible  only  for  a  discrete  set  of  values  of  k,  called  the  eigenvalues 
of  (7.4)  and  that  solutions  corresponding  to  different  eigenvalues  are 
orthogonal  to  each  other  in  the  sense 

/  *n  dv  =  0,  (7.6) 

whenever  kn  $  km.  Following  standard  practices  in  mathematical  physics, 
we  obtain  an  ascending,  infinite  sequence  of  eigenvalues 

kl  I  k2  1  k3  i  '• **  (7.7) 

and  a  corresponding  sequence  of  orthogonal  and  dimensionless  eigenfunctions 

^7  *2»  •  •  • 

so  normalized  that  they  satisfy  the  relations 

/  ^  <l»j  dv  =  5.^  V,  (7.8) 

where  dv  =  dx  dy  dz  and  integration  is  over  the  entire  volume  V  of  region 
ft.  It  is  known  that  the  sequence  of  such  functions  is  complete  for  equa¬ 
tions  of  the  type  we  a. a  considering  and  that  a  well-behaved  function  of 
x,  y,  z  is  expandable  in  the  region  ft  in  terms  of  these  eigenfunctions. 
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convergence  in  the  mean  is  assured  by  the  standard  theorems  of  functional 
analysis. 

With  these  preliminaries  we  can  reduce  the  solution  of  our  original 
problem  to  the  calculation  of  the  excitation  of  an  infinite  number  of 


oscillators  by  a  given  exciting  force.  We  write 

4  =  1  Un(t)<l>n(x,y,z),  (7.9) 

g  =  I  Kn(t)*n(x,y,z).  (7.10) 

The  coefficients  un(t)  and  g^(t)  are  calculated  as  follows: 

v  Wn  =  /  ♦♦n*  dv»  (7.11) 

V  Sn  =  /  g*n*  dv.  (7.12) 

The  differential  equation  (7.1)  then  requires  that 

un"  +  “n2  Wn  =  «n  (7.13) 


for  n  =  1,2,3,  *••.  Since  the  gn's  are  calculable  from  (7.12),  the  solu¬ 
tion  of  the  complete  problem  now  reduces  to  the  determination  of  the  ex¬ 
citations  of  the  harmonic  oscillators  whose  equations  are  given  above. 

The  excitations  depend  on  the  characteristic  frequencies  of  these  oscillators, 
wn,  and  on  the  functions  gj,(t). 

The  total  energy  fed  into  the  n-th  oscillator  is  proportional  to 

m 

/  l«n^t>l2<lt: 

o 

and  to  another  factor  which  represents  the  response  of  the  n-th  oscillator 
to  a  normalized  exciting  xorce.  This  response  factor  depends  on  uin  and  on 
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the  time  dependence  of  the  exciting  force  g^.  (It  is  to  be  noted  that  gR 

has  the  same  dimension  as  g,  which  is  not  that  of  force.  At  the  same  time, 
the  dimension  of  pn  is  that  of  $.) 

In  principle  the  calculations  require  the  solution  of  infinitely  many 
differential  equations  of  the  form  (7.13)  for  very  general  functions  gn(t) 

ot»  their  right  side.  Practically  the  situation  is  much  simplified  because 
the  physically  relevant  part  of  the  solution  may  be  obtained  by  solving 
only  one  or  a  fev  equations  of  the  type  described.  There  are  two  main 
reasons  for  this  simplification.  First,  in  a  situation  with  some  degree 
of  symmetry  many  of  the  functions  gn(t)  are  identically  zero,  because  the 
spatial  part  of  the  excitation  function  is  orthogonal  to  many  ij^'s.  Second, 
the  sequence  k^,  k2,  ***  increases  rather  rapidly  and  consequently  only  a 
few  of  the  harmonic  oscillators  which  are  described  by  equation  (7.13)  have 
resonant  frequencies  in  the  region  in  which  f  has  a  significant  harmonic 
component . 

He  shall  return  to  the  theory  after  exhibiting  some  of  the  well-known 
eigenvalues  and  eigenfunctions  appropriate  for  sound  propagation  in  a  long 
circular  cylinder. 

The  eigenfunctions  appropriate  for  a  right  circular  cylinder  of 
length  L  and  radius  b  are  the  non-vanishing  solutions  of  equation  (7.5) 
which  satisfy  the  boundary  conditions 

=  0  for  r  =  b,  and  =  0  for  z  =  ♦  L/2. 

3r  3z  — 

After  separating  equation  (7.5)  in  cylindrical  coordinates,  one  obtains 
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solutions  of  the  form 


$  =  e^n®  J  ( A‘2-h2r)  cos  hz  (7.14) 

s  n 

and 

=  ein6  Jn( A2-h2r)  sin  hz,  (7.15) 

where  the  fubscript  s  denotes  a  function  symmetric  for  the  reversal 
of  the  z-axis,  while  subscript  a  denotes  an  antisymmetric  function.  is 
the  Bessel  function  of  order  n,  where  n  must  be  an  integer  to  insure  unique¬ 
ness  of  the  solution.  These  functions  are  not  yet  normalized.  We  note  that 
if  the  cylinder  contained  one  or  more  septum;  i.e.  a  boundary  plane  contain¬ 
ing  the  cylinder  axis,  then  n  need  not  be  an  integer.  We  shall  avoid  this 
situation. 

The  possible  values  of  h  are  determined  by  the  boundary  conditions 
at  the  cylinder  ends.  These  require  that,  for  symmetric  functions  sin 
hL/2  =  0;  i.e. 

ht  =  2lir/L;  i  =  0,1,2,....  (7.16) 

For  antisymmetric  functions  cos  hL/2  =  0;  i.e. 

ht  =  (21  +  1)*/L;  t  =  0,1,2,....  (7.17) 

Boundary  conditions  on  the  cylinder  s' .  *  ,  require  that 

Jn'  ( A*-h2b)  =  0.  (7.18) 

Therefore  k  can  only  have  such  values  for  which  the  number  x  =  A^'-h^b 
is  a  solution  of  the  equation  Jn'(x)  =  0.  For  each  index  n,  the  deriva¬ 
tive  of  the  n-th  order  Bessel  function  has  an  infinite  number  of  discrete 

positive  roots.  We  write  these  in  ascending  order  as  x  ;  m  =  1,2,***. 

nm 

The  eigenvilues  of  the  differentia  .  equation  (7.5)  consistent  with  the 
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. . * . .  . . . . 


boundary  conditions  described  above  are  calculated  from  the  equations 

(kLl  “  hi)b2  =  xnm’  n  =  0,1,...;  m  =  1,2,...,  (7.19) 

where  h^  are  the  numbers  determined  from  (7.16)  and  (7.17).  The  symmetric 
eigenfunctions  are 

*snm£  8  ein9  Jn(xnmr/b)  cos  h£2’  (7*20> 

while  the  antisymmetric  eigenfunctions  are 

♦awnt  =  ftin6  Jn(xnmr/b)  sin  h£2*  (7*21) 

Of  the  triple  infinity  of  eigenfunctions  (or  oscillation  modes)  given  in 
(7.20)  and  (7.21),  certain  types  are  of  particular  interest  to  us.  Only 
£  =  0  produces  z-independent  eigenfunctions.  These  are  called  transverse 
modes.  Their  equations  are 

s  ein0  J  (X  r/b>.  (7.22) 

snm  n  nm 

Modes  which  have  rotational  symmetry  about  the  z-axis  are  all  of  the  form 

*som£  8  Jo(xomr/b)  cos  ^  <7-:3> 

and 

•  Jo‘W/b>  si"  (7-J4> 

The  first  few  roots  of  JQ'(r)  =  0  are  3.832,  7.016,  10.173,  13.324. 

Those  of  Ji'(r)  =  0  are  1.841,  5.330,  etc.11 

Let  us  now  calculate  the  excitation  of  a  cylindrical  resonator  with 
a  force  function  that  depends  only  on  r  and  t;  i.e.  which  is  independent 

II.  5ee  Ref.  8. 


-  38  - 


of  9  and  z.  In  this  case 


every  integral  of  the  type 


b  2n 

/  /  g(r,t 


o  o 


)  rdrd6 


vanishes  for  n  *  0,  leaving  only  the  eigenfunctions  with  n  =  0  to  interact 
with  g.  Similarly  the  integrals  over  z  involving  the  product  of  g  with 
sin  h£z  and  cos  h^z  will  vanish  except  for  cos  Oz.  Thus,  when  evaluating 
the  expansion  of  g(r,t)  in  terms  of  the  eigenfunctions  of  the  problem, 
only  the  symmetric  eigenfunctions  with  n  =  0  and  l  =  0  need  to  be  consid¬ 
ered.  The  relevant  eigenfunctions  are  J0(xomr/b).  Their  normalization 

factor  is  calculated  from  the  standard  formula12 
b 

l  CJ0(X0mr/b)]2rdr  =  A  b2[Jo(xom>J2,  (7.25) 


valid  for  every  xom  for  which  J0’(Xom)  =  0.  The  normalized  functions  are 

*m  =  Jo(W/b)/Jo<Xom).  (7.26) 

When  the  excitation  is  of  unit  magnitude  in  a  small  cylinder  of  radius  a 
centered  around  the  polar  axis  and  lasts  from  t  =  0  to  t  =  t,  we  have  from 
(7.12)  and  (7.26) 


Vgm(t)  =  J„(XL  )  f  J  <Xomr/b)rdrf(t),  (7.27) 

o '  *otn  o 

where  f(l)  =  1  sec"1  for  0  <  t  <  t,  and  0  for  t  >  t.  For  small  values  of 
a  we  nay  replace  J^X^r/b)  by  1  and  obtain  the  approximation 


*m(t) 


v  f(t) 

^  ^o^Xom)  * 


(7.28) 


where  v  =  a2irL  is  the  volume  swept  out  by  the  laser  beam. 


12.  Ref.  8,  p.9.l. 
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The  next  problem  is  to  deal  with  a  set  of  harmonic  oscillators  with 
pulse  excitation.  Before  losses  in  the  oscillators  are  considered,  their 
equations  have  the  form 

+  “2ym  =  *. .»  (7 .29) 

where  ^  is  given  in  (7.28).  The  characteristic  frequencies  are 

“m  *  X0mc/b-  (7.30) 

Taking  c  =  3.3  x  104  cm/sec,  we  get  for  the  lowest  symmetric  transverse 
frequency  (corresponding  to  xol  =  3.832) 

„  .  2.0xl04 

1  b 

and  for  the  next  lowest 

„  _  3.66xl04  , 

2  b 

wber®  b  cm»  v  i°  sec  Thus  the  lowest  such  frequency  for  a  tube 

10  cm  in  diameter  is  4000  sec*1,  the  next  7330  sec*1.  A  small  tube  of  1 
cm  diameter  would  require  acoustic  frequencies  which  are  prohibitively 
high. 

Generally  the  oscillators  will  have  a  certain  loss  rate  associated 
with  them.  Losses  are  not  only  unavoidable  but  necessary  when  the  res¬ 
onator  is  a  part  of  a  system  used  for  measurement,  because  measurement 
consists  of  conver* Ing  a  part  of  the  acoustic  disturbance  into  an  elec¬ 
tric  signal.  Thus  a  resonator  will  have  an  output  which  constitutes  the 
measured  signal,  or  else  the  energy  is  communicated  to  another  frequency- 
sensitive  device,  the  detector,  which  converts  it  to  an  electric  signal. 


In  any  case,  we  are  interested  in  devices  which  subject  the  disturbance 
to  a  harmonic  analysis.  The  noise  discriminating  property  of  the  detecting 
system  improves  as  its  bandwidth  is  restricted.  Unfortunately,  as  we  shall 
see,  narrowing  the  bandwidth  drastically  reduces  the  signal. 

For  the  sake  of  simplicity,  let  us  assume  that  the  detector  is  the 
cavity  itself,  that  is,  we  do  away  mentally  with  a  second  tuned  device  to 
which  the  cavity  may  be  coupled.  If  there  is  such  a  separate  device,  the 
detected  signal  will  2  less  because  some  of  the  energy  will  not  be  trans¬ 
ferred  and  new  sources  of  noise  will  be  introduced .  We  thus  have  a  sequence 
of  oscillators  with  an  output  and  we  concentrate  on  the  one  with  the  lowest 
frequency.  When  the  effects  of  oscillator  loss  (energy  conversion)  are  in¬ 
corporated  in  the  equations  of  the  harmonic  oscillators,  we  obtain  for  the 
oscillator  of  the  lowest  frequency13 

u"  +  Au'  +  (u^u  =  gi(t).  (7.31) 

The  quantity  A  has  the  dimension  of  reciprocal  time.  Let  W  denote  the 
energy  stored  in  the  free-running  harmonic  oscillator.  It  is  known  that 

A  =  w'-1  (7.32) 

dt 

One  can  readily  show  the  following:  When  g  is  a  sinusoidal  function 
of  time,  the  excitation  of  the  oscillator  is  significant  c.ily  if  the  fre¬ 
quency  of  the  excitirig  force  is  nearly  equal  to  that  of  the  free-runr.ing 
oscillator.  Then  by  means  of  Fourier  analysis  one  concludes  that  in  the 
case  of  a  general  time-dependent  force,  the  excitation  of  the  oscillator 

13.  P.M.  Morse,  Vibration  and  Sound,  McGraw-Hill,  New  York,  1948,  Ch.  II. 
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depends  largely  on  the  harmonic  component  of  the  exciting  .orce  which  lies 


near  v^.  After  a  short  calculation  one  finds  tnat  the  average  energy  of 

the  oscillator  is  approximately 

12 


W  % 


2A 


(7.33) 


where  G(v)  is  the  Fourier  transform  of  the  exciting  function  g(t)  and 
V1  "  wi/2*»  The  power  output  of  such  an  oscillator  is 


dW 

dt 


=  AW  i-lGtvi)!2 


(7.34) 


Only  part  of  this  is  useful  power  output  in  the  detector,  but  this  brings 
in  only  a  factor  of  about  one-half,  which  is  not  significant.  The  important 
matter  is  that  the  rate  at  which  power  is  absorbed  by  the  oscillators  is 
proportional  to  the  Fourier  transform  of  the  excitation  function  at  the 
resonant  frequencies  of  the  oscillators. 


The  effect  of  the  confining  tu^e  is  that  the  acoustic  signal  generated 
within  it  is  resolved  into  its  spectral  components  according  to  the  resonant 
frequencies  of  the  cavity.  As  we  have  seen  in  section  6,  the  power  spec¬ 
trum  of  a  single  pulse  of  duration  t  io  rather  broad,  most  of  the  energy  is 
in  the  spectral  interval  0  to  1/t.  Periodic  repetition  of  the  pulse  pro¬ 
duce-  a  spectral  distribution  of  the  acoustic  energy  among  the  harmonics  of 
the  repetition  frequency  in  the  manner  illustrated  in  rig.  3. 


In  o.xler  to  obtain  good  utilization  of  the  acoustic  energy  generated, 
the  experiment  should  be  so  designed  that  a  large  fraction  cf  the  acoustic 
energy  generated  by  tie  source  is  concentrated  in  the  frequency  region  cap¬ 
able  of  exciting  a  single  selected  cavity  resonance  because  the  detector 
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will  only  be  coupled  to  one  resonance.  Such  concentration  of  the  acoustic 
energy  is  not  easily  achieved.  With  a  single  pulse  it  is  just  about  hope¬ 
less  to  obtain  success  in  this  manner.  With  pulses  repeated  at  the  rate  of 
a  resonant  frequency  of  the  cavity,  the  situation  is  better,  but  reason¬ 
able  cavity  dimensio  s  require  rather  high  pulse  repetition  frequencies. 

We  calculated  that  the  lowest  resonance  of  a  cylindrical  cavity  of  10  cm 
diameter  is  at  4000  Hz,  therefore  we  would  have  to  operate  with  the  laser 
pulses  shorter  than  0.1  msec  repeated  4000  times  per  second.  The  success 
of  the  scheme  of  this  type  depends  not  only  on  the  construction  of  a  power¬ 
ful  laser  with  appropriate  pulsing  mechanism  but  also  on  :he  construction 
of  a  highly  sensitive  acoustic  detector  coupled  with  the  proper  mode  of 
oscillation  of  a  cylinder. 

8.  Conclusions. 

The  following  conclusions  were  reached  concerning  the  measurement  of 
the  absorption  cf  radiation  in  air  when  tb-  wavelength  of  the  radiation  is 
in  the  0.5  to  1.2  micron  range: 

1.  It  is  probably  possible  to  make  absorption  measurements  in  air 
using  a  steady  laser  beam  an'4  noting  the  static  pressure  in¬ 
crease  in  a  sealed  tube  if  the  following  conditions  can  be  met: 

a)  The  tube  is  long  enough  to  render  the  end  effects  negligible. 

The  length  should  be  at  least  5000  times  the  thickness  of  a 
window. 

b)  The  heating  effect  of  the  scattered  light  is  eliminated. 

c)  The  temperature  of  the  confining  tube  is  maintained  constant 
to  within  10"4  °C. 
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2.  The  static  pressure  rise  due  to  a  single  laser  pulse  can  probably 
not  be  measured  if  a  is  as  low  as  10-9  cm-1. 

3.  The  measurement  of  the  acoustic  wave  generated  by  a  single  laser 
pulse  depends  on  the  focussing  of  the  acoustic  wave  and  on  an 
extremely  sensitive  broad-band  detector.  It  is  doubtful  that 
this  can  be  achieved. 

4.  The  measurement  of  the  acoustic  signal  generated  by  periodically 
repeated  pulses,  whose  repetition  frequency  is  matched  to  that 

of  a  resonant  cavity,  is  a  possibility  that  should  be  investigated 
further. 

The  situation  in  the  wavelength  region  around  10  microns  is  quite  different 
from  that  investigated  here.  First,  Rayleigh  scattering  is  no  longer  a 
limiting  factor.  Second,  the  absorption  in  air  might  be  considerably  higher 
at  10  microns  than  it  is  in  the  visible  region.  If  this  should  be  the  case, 
several  of  the  contemplated  methods  would  work. 

~iiis  work  does,  not  encomoass  the  evaluation  of  optical  metnods  which 
should  i»e  studied  separately. 
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Figure  1 


.  Distortion  of  a  square  pulse  emanating  'om  a  line  source. 
The  function  ♦(R,t)  defined  by  equation  (6.40)  is  shown 
for  two  different  values  of  R.  (R'  =  R/c.) 


Figure  2.  The  spectral  density  function  of  a  square  pulse. 
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require,  provided  that  in  either  case  the  release,  use, 
duplication  or  disclosure  hereof  shall  be  subject  to  the 
foregoing  limitations.  This  legend  shall  be  marked  on 
any  reproduction  hereof  in  whole  or  in  part." 


If  the  above  statement  is  to  be  used  on  this  form  enter 
the  following  abbreviated  statement: 

"Furnished  under  U.  S.  Government  Contract  No. _ . 

S  tall  not  be  either  released  outside  the  Government,  or  used, 
duplicated,  or  disclosed  in  whole  or  in  part  for  manufacture 

or  procurement,  without  the  written  permission  of _ . 

per  ASPR  9-203.  " 

DoD  Imposed  D  stnbution  Statement  reference  DoD  /Jire.  five 
5 200  2o  )  "Distribution  Statements  ( other  than  Security)  on 
Technical  Documents,"  March  29.  1965. 

STATEMENT  NO.  1  -  Distribution  of  this  document  is  unlimited. 

STATEMENT  NO  2  (UNCLASSIFIED  document)  -  This  document 
ia  subject  to.  special  export  controls  snd  each  transmittal  to 
foreign  governments  or  foreign  nationals  may  be  made  only  with 
prior  approval  of  (/ill  in  controlling  DoD  office). 

(CLASSIFIED  document)  •  In  add. non  to  aeevrity  require¬ 
ments  which  must  be  met,  this  document  ia  subject  to  special 
export  controls  and  each  transmittal  to  foreign  governments  or 
foreign  nationals  may  be  made  only  with  prior  approval  (fill 
in  controlling  DoD  Office). 

STATEMENT  NO  3  (UNCLASSIFIED  document)  -  Each  trans¬ 
mittal  of  this  document  outside  the  agencies  of  the  U.  S. 
Government  must  have  prior  approval  of  (fill  in  controlling 
DoD  Office). 

(CLASSIFIED  document)  -  In  addition  to  security  require¬ 
ments  which  apply  to  thia  document  and  must  be  met,  each 
transmittal  eutaide  the  agenclea  of  the  U.  S.  Government  must 
have  prior  approval  of  (fill  in  controlling  DoD  Office). 

STATEMENT  NO  4  (UNCLASSIFIED  document)  -  Each  trans¬ 
mittal  of  this  document  outside  the  Department  of  Defense 
must  have  prior  approval  of  (HI!  in  controlling  DoD  Office ). 

(CLASSIFIED  document)  -  In  addition  to  security  require¬ 
ments  which  apply  to  this  document  and  must  be  met,  each  trana- 
mittsl  outside  the  Department  of  Defense  must  have  prior  ap¬ 
proval  of  (fill  i' i  controlling  DoD  Office). 

STATEMENT  NO.  S  (UNCLASSIFIED  document)  -  This  document 
may  be  further  distributed  by  cny  holder  only  with  specific 
prior  approval  of  (till  m  controlling  DoD  Office). 

(CLASSIFIED  document)  •  In  addition  to  security  require¬ 
ments  which  apply  to  this  document  and  must  be  met.  it  may 
Ire  further  distributed  b>  the  holder  ONLY  with  specific  prior 
approval  of  (I'll  in  controlling  DoD  Oftu  e). 

11.  SUPPLEMENTARY  NOTES.  Use  for  additional  explanatory 
notea. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponaoring  (paying 
for)  the  research  and  development.  Include  addreas. 

12.  ABSTRACT:  Enter  an  abstract  givinga  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  t1  body  of  the  technical  re¬ 
port.  If  additional  space  is  req"  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  Ihc  abstract  of  classified  re¬ 
porta  be  unclassified.  Each  paragraph  of  the  abstract  shall 
end  with  an  indication  of  the  military  security  classification 
of  the  information  in  the  paragraph,  represented  as  (75).  fS). 

(C).  or  {(/> 


There  ia  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Kay  words  must  be 
selected  so  that  no  secutity  claaaification  is  required.  Iden¬ 
tifiers,  such  as  equipment  model  designation,  trad-r  name, 
military  project  code  name,  geographic  location.  May  be  u:,ed 
sa  key  words  but  will  be  followed  by  an  indication  of  techni¬ 
cal  context.  The  assignment  of  links,  roles,  and  weights  ia 

optional.  D*rc  it 
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